Tunneling currents in ferromagnetic systems with multiple broken symmetries 
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A system exhibiting multiple simultaneously broken symmetries offers the opportunity to influence physical 
phenomena such as tunneling currents by means of external control parameters. In this paper, we consider the 
broken SU (2) (internal spin) symmetry of ferromagnetic systems coexisting with i) the broken (7(1) symmetry 
of superconductors and ii) the broken spatial inversion symmetry induced by a Rashba term in a spin-orbit cou- 
pling Hamiltonian. In order to study the effect of these broken symmetries, we consider tunneling currents that 
arise in two different systems; tunneling junctions consisting of non-unitary spin-triplet ferromagnetic super- 
conductors and junctions consisting of ferromagnets with spin-orbit coupling. In the former case, we consider 
different pairing symmetries in a model where ferromagnetism and superconductivity coexist uniformly. An 
interplay between the relative magnetization orientation on each side of the junction and the superconducting 
phase difference is found, similarly to that found in earlier studies on spin-singlet superconductivity coexisting 
with spiral magnetism. This interplay gives rise to persistent spin- and charge-currents in the absence of an 
electrostatic voltage that can be controlled by adjusting the relative magnetization orientation on each side of 
the junction. In the second system, we study transport of spin in a system consisting of two ferromagnets with 
spin-orbit coupling separated by an insulating tunneling junction. A persistent spin-current across the junction is 
found, which can be controlled in a well-defined manner by external magnetic and electric fields. The behavior 
of this spin-current for important geometries and limits is studied. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.20.-z 



I. INTRODUCTION 

Due to the increasing interest in the field of spintronics in 
recent years 1 , the idea of utilizing the spin degree of free- 
dom in electronic devices has triggered an extensive response 
in many scientific communities. The spin-Hall effect is ar- 
guably the research area which has received most focus in 
this context, with substantial effort being put into theoreti- 
cal considerations 2 as well as experimental observations 3 -. In 
spintronics, a main goal is to make use of the spin degree of 
freedom rather than electrical charge, investigations of mech- 
anisms that offer ways of controlling spin-currents are of great 
interest. The study of systems with multiple broken sym- 
metries is highly relevant in this context, since such systems 
promise rich physics with the opportunity to learn if the tun- 
neling currents can be influenced by means of external control 
parameters such as electric and/or magnetic fields. Here, we 
will focus on two specific systems: ferromagnetism coexisting 
with superconductivity, which we shall refer to as ferromag- 
netic superconductors (FMSC), and systems where ferromag- 
netism and spin-orbit coupling are present (FMSO). In terms 
of broken symmetries, we will then study the broken SU (2) 
(internal spin) symmetry of ferromagnetic systems coexisting 
with the broken U(l) symmetry of superconductors and also 
consider ferromagnets with broken inversion (spatial) symme- 
try induced by a Rashba term in a spin-orbit coupling Hamil- 
tonian. 

The coexistence of ferromagnetism (FM) and superconduc- 
tivity (SC) has a short history in experimental physics 4,5,6 , 
although a theoretical proposition of this phenomenon was 
offered as early as 1957 by Ginzburg 7 . Spin-singlet super- 
conductivity originating with BCS theory seems to be ruled 
out as a plausible pairing mechanism for a ferromagnetic 



superconductor , at least with regard to uniform coexistence 
of the FM and SC order parameters £ and A, respectively. It 
could be achieved for a superconductor taking up a so-called 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 9 . However, 
it seems likely that the coexistence of FM and SC call foriSJi 
p-wave spin-triplet Cooper pairs which have a non-zero mag- 
netic moment. This type of pairing has been observed in 
superfluid 3 He, and is perfectly compatible with FM order. 
Spin-triplet superconductivity has moreover been experimen- 
tally verifie d 1213 in Si"2RuC»4, and the study of such a pair- 
ing in a FMSC could unveil interesting effects with respect 
to quantum transport. The concept of simultaneously broken 
U(l) and SU(2) symmetries are of great interest from a fun- 
damental physics point of view, and could be suggestive to a 
range of novel applications. This topic has been the subject of 
theoretical research in e.g. Refs. 14. 1 5lfl6l 

In this paper, we follow up Ref. [Tt| with a more compre- 
hensive study of the tunneling currents between two p-wave 
FMSC separated by an insulating junction; RuS^GdQ^Os, 
UGe2, and URhGe have been proposed as candidates for such 
unconventional superconductors^ 6 -. In our model, we as- 
sume uniform coexistence of the FM and SC order parame- 
ters and that superconductivity arises from the same electrons 
that are responsible for the magnetism. As argued in Ref. HI 
this can be understood most naturally as a spin-triplet rather 
than spin-singlet pairing phenomenon. Furthermore, it seems 
that SC in the metallic compounds mentioned above always 
coexists with the FM order and is enhanced by ifcl£; the ex- 
periments conducted on the compounds UGe2 and URhGe do 
not give any evidence for the existence of a standard normal- 
to-superconducting phase transition in a zero external mag- 
netic field, but instead indicate a phase corresponding to a 
mixed state of FM and SC. We provide detailed calculations 
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for single-particle and Josephson (two-particle) tunneling be- 
tween two non-unitary equal-spin pairing (ESP) FMSC. We 
examine both the charge- and spin-sector in detail within lin- 
ear response theory using the Kubo formula. We find that the 
supercurrent of spin and charge may be controlled by adjust- 
ing the misorientation of the exchange fields on both sides of 
the junction. Such an effect was first discovered by Kulic and 
Kulic 19 , who derived an expression for the Josephson current 
over a junction separating two BCS superconductors with spi- 
ral magnetic order. It was found that the supercurrent could be 
controlled by adjusting the relative orientation of the exchange 
field on both sides of the junction, a finding that quite re- 
markably suggested a way of tuning a supercurrent in a well- 
defined manner from e.g. a 0- to ir -junction. Later investiga- 
tions made by Eremin, Nogueira, and Tarentc 20 considered 
a similar system as Kulic and Kulic 19 , namely two Fulde- 
Ferrel-Larkin-Ovchinnikov (FFLO) superconductors 9 coex- 
isting with helimagnetic order. Recently, the same opportu- 
nity was found to exist in a FMSC/I/FMSC junction as shown 
by Gr0nsleth et alM. 

In the case of a system where both ferromagnetism and 
spin-orbit coupling are present, it is clear that these are physi- 
cal properties of a system that crucially influence the behavior 
of spins present in that system. For instance, the presence 
of spin-orbit coupling is highly important when considering 
ferromagnetic semiconductors 21,22 . Such materials have been 
proposed as devices for obtaining controllable spin injection 
and manipulating single electron spins by means of external 
electrical fields, making them a central topic of semiconduc- 
tor spintronics 2 ™ In ferromagnetic metals, spin-orbit coupling 
is ordinarily significantly smaller than for semiconductors due 
to the bandstructure. However, the presence of a spin-orbit 
coupling in ferromagnets could lead to new effects in terms of 
quantum transport. 

Studies of tunneling between ferromagnets have uncovered 
interesting physical effect s 24 ' 25 ' 26 . Nogueira et al. predicted 24 
that a dissipationless spin-current should be established across 
the junction of two Heisenberg ferromagnets, and that the 
spin-current was maximal in the special case of tunneling be- 
tween planar ferromagnets. Also, there has been investiga- 
tions of what kind of impact spin-orbit coupling constitutes 
on tunneling currents in various contexts, e.g. for noncen- 
trosymmetric superconductors 27 , and two-dimensional elec- 
tron gases coupled to ferromagnets 28 . Broken time reversal- 
and inversion-symmetry are interesting properties of a system 
with regard to quantum transport of spin and charge, and the 
exploitment of such asymmetries has given rise to several de- 
vices in recent years. For instance, the broken SU(2) sym- 
metry exhibited by ferromagnets has a broad range of possi- 
ble applications. This has led to spin current induced mag- 
netization switching 29 , and suggestions have been made for 
more exotic devices such as spin-torque transistors 30 and spin- 
batteries 31 . It has also led to investigations into such phe- 
nomena as spin-Hall effect in paramagnetic metals 32 , spin- 
pumping from ferromagnets into metals, enhanced damp- 
ing of spins when spins are pumped from one ferromagnet 
to another through a metallic sample^ 3 ., and the mentioned 
spin Josephson effects in ferromagnet/ferromagnet tunneling 



junctions* 4 -. 

Here, we study the spin-current that arises over a tunneling 
junction separating two ferromagnetic metals with substantial 
spin-orbit coupling. It is found that the total current consists 
of three terms; one due to a twist in magnetization across the 
junction (in agreement with the result of Ref. [24]), one term 
originating from the spin-orbit interactions in the system, and 
finally an interesting mixed term that stems from an interplay 
between the ferromagnetism and spin-orbit coupling. After 
deriving the expression for the spin-current between Heisen- 
berg ferromagnets with substantial spin-orbit coupling, we 
consider important tunneling geometries and physical limits 
of our generally valid results. Finally, we make suggestions 
concerning the detection of the predicted spin-current. Our 
results indicate how spin transport between systems exhibit- 
ing both magnetism and spin-orbit coupling can be controlled 
by external fields, and should therefore be of considerable in- 
terest in terms of spintronics. 

This paper is organized as follows. In Sec. [II] we consider 
transport between spin-triplet ferromagnetic superconductors, 
while a study of transport between ferromagnets with spin- 
orbit coupling is given in Sec. [Ill] A discussion of our results 
is provided in Sec. IIVI with emphasis on how the novel ef- 
fects predicted in this paper could be tested in an experimental 
setup. Finally, we give concluding remarks in Sec. [V] 



II. FERROMAGNETIC SUPERCONDUCTORS 

A. Coexistence of ferromagnetism and superconductivity 

An important issue to address concerning FMSC is whether 
the SC and FM order parameters coexist uniformly or if 
they are phase-separated. One possibility 34 is that a sponta- 
neously formed vortex lattice due to the internal magnetiza- 
tion m is realized in a spin-triplet FMSC, while there also 
have been studies of Meissner (uniform) SC phases in spin- 
triplet FMSGi£. As argued in Ref.[35l a key variable with re- 
spect to whether a vortice lattice appears or not is the strength 
of the internal magnetization m. Ref. [36] suggested that vor- 
tices arise if 47rm > H c i, where H c i is the lower critical 
field. When considering a weak FM state coexisting with SC, 
a scenario which seems to be the case for URhGe, the domain 
structure in the absence of an external field is thus vortex-free. 
Current experimental data concerning URhGe are not strong 
enough to unambiguously settle this question, while evidence 
for uniform coexistence of FM and SC has been indicated 37 
in UGe2. Furthermore, a bulk Meissner state in the FMSC 
RuSr^GdCu^Os has been reported in Ref.[38l hence suggest- 
ing the existence of uniform FM and SC as a bulk effect. In 
our study, we shall consequently take the order parameters as 
coexisting homogeneously and use their bulk values, as justi- 
fied by the argumentation above. However, we emphasize that 
one in general should take into account the possible suppres- 
sion of the SC order parameter in the vicinity of the tunnel- 
ing interface due to the formation of midgap surface states 39 
which occur for certain orientations of the SC gap. The pair- 
breaking effect of these states in unconventional superconduc- 
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tors has been studied in e.g i 40,41,42 , and we discuss this in more 
detail in Sec. [IV] A sizeable formation of such states would 
suppress the Josephson current, although it is nonvanishing 
in the general case. Also, we use bulk uniform magnetic or- 
der parameters, as in 24 . The latter is justified on the grounds 
that a ferromagnet with a planar order parameter is mathemat- 
ically isomorphic to an s-wave superconductor, where the use 
of bulk values for the order parameter right up to the interface 
is a good approximation due to the lack of midgap surface 
states. 

It is generally believed that the same electrons that are re- 
sponsible for itinerant FM also participate in the formation 
of Cooper pairs below the SC critical temperature 6 . As a 
consequence, uniform coexistence of spin-singlet SC and FM 
can be discarded since s-wave Cooper pairs carry a total spin 
of zero, although spatially modulated order parameters could 
allow for magnetic s-wave superconductors 1 ^. However, 
spin-triplet Cooper pairs are in principle perfectly compati- 
ble with FM order since they can carry a net magnetic mo- 
ment. To see this, consider the d k -vector formalism 43 which 
is convenient when dealing with spin-triplet superconductors, 
regardless of whether they are magnetic or not. For a complete 
and rigorous treatment of the dk-vector order parameter, see 
e.g. Ref. |44j. The spin dependence of triplet pairing can be 
represented by a 2 x2 matrix 



AkiT AkU. 
—d x (k) +id y (k) d z 
d z (k) d x (k)- 



(1) 



(k) 

Md y (k) 



= ld k • (T(T y 



where A kQ/ 3 represent the SC gap parameters for different 
triplet pairings, a = (a x , a y ,a z ) where &i are the Pauli ma- 
trices, and dk = (d x (k), d y (k), d z (k)) is given by 



/ Akxi - A kTT 
V 9. 



Men 



A, 



Note that dk transforms like a vector under spin rotations and 
that Akjx = Ak|| for triplet pairing since it is of no signif- 
icance which electron in the Cooper pair that has spin up or 
down. This is because spin-part of the two-particle wavefunc- 
tion is symmetric under exchange of particles, as opposed to 
spin-singlet SC, where the gap changes sign when the spin in- 
dices are exchanged. Spin-triplet SC states are classified as 
unitary if idk x d k = and non-unitary if the equality sign 
does not hold. Since the average spin of a dk-state is given 



(Sk) 



id k x d k , 



(3) 



it is clear that we must have a non-unitary dk in a model where 
FM and SC coexist uniformly. Indeed, there is strong reason 
to believe that the correct pairing symmetries in the discovered 
FMSC constitute non-unitary states 1 'i 45 ? 46 . As a consequence, 
one can rule out for instance a state where only Akfj, ^ 
since it would imply (Sk) = according to Eq. (0. In the 
most general case where all SC gaps are included, Ajx would 
be suppressed in the presence of a Zeeman-splitting between 
the |, J. conduction bands 6 ; see Fig. [T] 




FIG. 1: Band-splitting for f, J. electrons in the presence of a magne- 
tization in z-direction. Inter-band pairing gives rise to a net Cooper 
pair momentum in the presence of a band-splitting, thus suppressing 
the AkTi order parameter. 



However, such a splitting between energy-bands need not 
be present and one could in theory then consider a d-vector 
where 



+ 0, A kU ± 



(4) 



such that (Sk) lies in the local xj/-plane. This scenario would 
be equivalent to an A2 -phase as is seen when performing a 
spin rotation on the gap parameters into a quantization axis 
lying in the xy-plane. Denoting up- and down-spins with re- 
spect to the new quantization axis by + and — , respectively, 
the transformation yields 




1 

e -i0 

-2nj5 



2e 14 


-2e~ 




, (5) 



kTJ-J- ^ where is the azimuthal angle as shown in Fig. |2] 



® 



*kU 




FIG. 2: Change of spin-basis for the superconducting gaps. The new 
quantization axis is represented by the dotted arrow. 

When introducing the conditions in Eq. ©, it is readily 
seen that A k+ _ = while |A k+ +| ^ |A k | ^ 0, thus 
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corresponding to an A2 -phase. Consequently, the entire span 
of physically possible pairing symmetries in a FMSC can be 
reduced to the equivalence of an Al- or A2-phase in 3 He by 
a change of spin-basis. The definitions of A-, A1-, and A2- 
phases in 3 He are as follow: an A-phase corresponds to a pair- 
ing symmetry such that |Ak|t — I Aim I 0, an Al-phase 
has only one gap A-^ aa ^ while A^-o-.-er = 0, and an 
A2-phase satisfies |Akf-f | ^ |Akxxl ^ 0. In this case, Ak Q( 3 
represents the superfiuid gap for the fermionic 3 He-atoms, and 
Ak|x = for all ^-phases. 

The resulting spin of the Cooper pair is then in general 
given by 



Hamiltonian similar to the one used in Ref . |4£ 



(Sk} = (l/2)[|A kTT | 2 -|Aku| 2 ]z. 



(6) 



In the following, we shall accordingly consider tunnel- 
ing between non-unitary ESP FMSC in an Al- or A2-phase. 
Moreover, we consider thin film FMSC, ensuring that no ac- 
cumulation of charge at the surface will take place due to an 
orbital effect. Our system can be thought to have arisen by 
first cooling down a sample below the Curie temperature Tm 
such that FM order is introduced. At further cooling below 
the critical temperature T c , the same electrons that give rise to 
FM condense into Cooper pairs with a net magnetic moment 
parallel to the original direction of magnetization. Our model 
is shown in Fig. [3] 




FIG. 3: Tunneling between two non-unitary ESP FMSC. The quan- 
tization axis has been taken along the direction of magnetization on 
each side of the junction. 



B. The Hamiltonian 

The system consists of two FMSC separated by an insulat- 
ing layer such that the total Hamiltonian can be written as 47 
H = Hl + iJp> + Hj, where L and R represents the indi- 
vidual FMSC on each side of the tunneling junction, and Hj 
describes tunneling of particles through the insulating layer 
separating the two pieces of bulk material. Using mean-field 
theory, one finds that the individual FMSC are described by a 



-Hfmsc = H + 2_^ ip^Akipk, 

k 

H = JN V (0)m 2 + X - J2 £k. + A U W ( ? > 

kcr ka/3 

Here, k is the electron momentum and we have introduced 

Ek CT =s k -aQ, o-=T,|=±l. (8) 

Furthermore, J is a spin coupling constant, 77 (k) is a geomet- 
rical structure factor which for k = reduces to the number 
of nearest lattice neighbors 77(0), m = {m x , m v , m z } is the 
magnetization vector, while Ak Q /3 is the superconducting or- 
der parameter and b-^ a p = (c_k/3Ck Q ) denotes the two-particle 
operator expectation value. The ferromagnetic order parame- 
ters are given by 



C = 2Ji](Q)(m x - vm y ), C, z = 2Jr)(0)m z 



(9) 



The interesting physics of the FMSC/FMSC junction lies in 
the matrix „4k to be given below. Above, we used a basis 



i'k = (ckt Cki c^ kT cL ki ) 



(10) 



where Ck CT (c ko .) are annihilation (creation) fermion operators. 
Note that we have not incorporated any spin-orbit coupling of 
the type (E x p) • cr in the Hamiltonian described in Eq. (0 
such that spatial inversion symmetry is not broken, i.e. we 
consider centrosymmetric FMSC. 
Consider now the matrix 



/-£k T C A 



■Ak 



ki r 



A v 



C f -Ekj. A kiT A kU 
-C e n ) 



A ku 



(ID 



which is valid for a FMSC with arbitrary magnetization. As 
explained in the previous sections, we will study in detail tun- 
neling between non-unitary ESP FMSC, i.e. Akix = Akxx = 
0, C = in Eq. dTTb . We take the quantization axis on each 
side of the junction to coincide with the magnetization direc- 
tion. One then needs to include the Wigner rf-function 4 ^ de- 
noted by fP) a (€) with j = 1/2 to account for the fact that 
a | spin on one side of the junction is not the same as a j 
spin on the other side of the junction, since the magnetization 
vectors point can point in different directions. The angle d is 
consequently defined by 



iiir • rrii = m^m^ cos(z9), m, t 
Specifically, we have that 



■n(i/a)w_ fcos^/2) -sin(tf/2)\ 
U [V) ~ ^sin(i?/2) cos(^/2) J 



(12) 



(13) 
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such that a spin-rotated fermion operator is given by 



(14) 



The tunneling Hamiltonian then reads 

kpcrcr' 

where we neglect the possibility of spin-flips in the tunneling 
process. Note that we distinguish between fermion operators 
on the right and left side of the junction corresponding to Cka 
and d pa , respectively. Demanding that Hj is invariant under 
time reversal /C, one finds that the condition Kr l HiK. = Hj 
with 



K~ 1 H T K = £ aa'V^Wx 

kpcrcr' 

( T kpcL k: _ CT cLp,-<T' H 



(16) 



dictates that Tkp = T*k p - Furthermore, we write the super- 
conducting order parameters as Ako-o- = |Al C o- C r|e 1 ^ k+9 ''< T '', 
where R (L) denotes the bulk superconducting phase on the 
right (left) side of the junction while 6*k is a general (complex) 
internal phase factor originating from the specific form of the 
gap in k-space that ensures odd symmetry under inversion of 
momentum, i.e. #k = #-k + 7r- 

For our system, Eq. (O takes the form 



FMSC 



H + H A , H 



(17) 



C. Tunneling formalism 

Although the treatment in this section is fairly standard, it 
comes with certain extension to the standard cases due to the 
coexistence of two simultaneously broken symmetries. Thus, 
for completeness, we present it here. 

In order to find the spin- and charge-current over the junc- 
tion, we define the generalized number operator 67 by N a p — 
Sk c ka c k/3- Consider now the transport operator 



N a p =i[Ht, N a/3 ] 

kpcr 



pcr C kp\ 



(22) 

We now write H = H' + H T where H' = H L + H R and 
Hi — Ki + /ijA/j, i = L, R, where fn is the chemical potential 
on side i and Ni is the number operator. In the interaction 
picture, the time-dependence of N a p is then governed by 



N af3 (t)=e lH ' t N a0 e- lH ' t , 



(23) 



while the time-dependence of the fermion operators reads 

CkaW = e lK ^c ka e- lK ^. (24) 
Effectively, one can write 

K R = H + E #k CT 7k CT 7k CT , (25) 



where we have block-diagonalized A k and chosen a conve- 
nient basis k(T = (c k(r , c_kcr), with the definition 



A ka = -- 



1 /-e kCT A 



kcrcr 



2 VAL. 



(18) 



This Hamiltonian is diagonalized by a 2 x 2 spin generalized 
unitary matrix U krJ , so that the superconducting sector is ex- 
pressed in the diagonal basis 



*L = ^kaUka = (7 k(7 >7- kCT )- 



(19) 



Thus, H A = Ska ^L^k CT 0ker, in which 



= Uk*A ka U k J: = diag(-E k(7 , -E ka )/2, 



E ^ = \I £ L + l A k<x<x| 2 - 
The explicit expression for U ka is 

1 



(20) 



iVka = 



Ekcr+Skcr 
£k(T + ^k(T 



(21) 



(£kcr + -Eko 



ik(T(T 



We now proceed to investigate the tunneling currents that can 
arise across a junction of two such FMSC. 



where the chemical potential is now included in the quasi- 
particle excitation energies E ka according to 



(26) 



with £kcr = £kcr — /-tR, and correspondingly for the left side. 
Consequently, we are able to write down 

N a p(t) = -iE (^ /2) W^ P cLW^(t)e- iteV 

k P * (27) 

-V^' 2 \d)T^dl a {t)c^ teV ), 

where eV = /Ul — A*R i s tr,e externally applied potential. 
Within linear response theory, we can identify a general cur- 
rent 

I(t) = E r a p{N^{t)), t = (-el, <x), (28) 

a/3 

such that the charge-current is 7 C (t) = Jo if) while the spin- 
current reads I s (t) = (Ix(t),I 2 (t),I 3 (t)). In Eq. ([28), 1 
denotes the 2x2 identity matrix. Explicitely, we have 

/ c (^j s c p W + /JW = - e ^(iv Qa (t)) 

a 

I s (i) = I%(t) +If p (t) =E^(^(*))= (29) 

a/3 
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where the subscripts sp and tp denote the single-particle and 
two-particle contribution to the currents, respectively. As 
recently pointed out by the authors of Ref. 50, defining a 
spin-current is not as straight-forward as defining a charge- 
current. Specifically, the conventional definition of a spin- 
current given as spin multiplied with velocity suffers from 
severe flaws in systems where spin is not a conserved quan- 
tity. In this paper, we define the spin-current across the junc- 
tion as I s (t) = (dS (*)/&) where dS/dt = i[# T ,S]. It is 
then clear that the concept of a spin-current in this context 
refers to the rate at which the spin-vector S on one side of the 
junction changes as a result of tunneling across the junction. 
The spatial components of I s are defined with respect to the 
corresponding quantization axis. In this way, we avoid non- 
physical interpretations of the spin-current in terms of real 
spin transport as we only calculate the contribution to dS /dt 
from the tunneling Hamiltonian instead of the entire Hamil- 
tonian H. Should we have chosen the latter approach, one 
would in general run the risk of obtaining a non-zero spin- 
current due to e.g. local spin-flip processes which are obvi- 
ously not relevant in terms of real spin transport across the 
junction. However, in our system such spin-flip processes are 
absent. 

The tunneling currents are calculated in the linear response 
regime by using the Kubo formula, 



(N af3 (t)) = -i / dt'([N a p(t),H T (t')}), 



(30) 



where the right hand side is the statistical expectation value in 
the unperturbed quantum state, i.e. when the two subsystems 
are not coupled. This expression includes both single-particle 
and two-particle contributions to the current. Details of the 
calculations are found in Sec. I A II 



We now consider the cases of an A2- and A 1 -phase at 
zero external potential, giving special attention to the charge- 
current and z-component of the spin-current in the Josephson 
channel. 



with the definitions 

I C (<t>u<t>2) = |5}l + <7<*cos(tf)]|r kp 



er^,. _ , 2 |A kaQ ||A pcrcr | 



x cos(^i) sm(<p 2 )F kpa(T , 
/ S (<^ 2 ) = -^£c*[l + ™ «»(*)] |T kp 1 2 |Ak ; al1 ^ 1 

x cos(0i) sin(0 2 )i ; kpa ( T, (32) 
where we have introduced A6> pk = P — 6* k and 



Fi 



kpacr 



E 



f(±E ka ) - f(E pa ) 



(33) 



Above, f(x) is the Fermi distribution. Thus, we have found 
a two-particle current, for both spin and charge, that can be 
tuned in a well-defined manner by adjusting the relative ori- 
entation -d of the magnetization vectors 68 . We will discuss the 
detection of such an effect later in this paper. Note that the k- 
dependent symmetry factor 6> k enters the above expressions, 
thus giving rise to an extra contribution to the two-particle 
current besides the ordinary Josephson effect. This is due to 
the fact that we included it in the SC gaps as a factor e k 
which in general is complex. However, this specific form may 
for certain models, depending on the Fermi surface in ques- 
tion, be reduced to a real function, i.e. e l6k — > cos(9 k , in 
which case the quasi-particle interference term becomes zero. 
Hence, in most of the remaining discussion we will focus on 
the Josephson part of the two-particle current. 

The A 1 -phase with only one SC order parameter A kQQ , 
a e {I, [} also corresponds to a non-unitary state d k accord- 
ing to Eq. (01, and is thus compatible with coexistence of FM 
and SC. In this case, we readily see that Eq. (f32t reduces to 



lf v = ecos 2 (d/2)X a 
-acos 2 (d/2)X 



6{T,1} 



(34) 



tp.z 



where we have defined the quantity 



D. Two-particle currents 

For an A2 -phase in the case of zero externally applied volt- 
age (eV = 0), Eqs. ( f29b and jAl 3I > generates a quasiparti- 
cle interference term I q \, in addition to a term Ij identified as 
the Josephson current. Thus, the total two-particle currents 
of charge and spin can be written as I^fl^ 
where 



jC(S) 

qi(>z) 



rC(S) 



C)=E /C(S) (^-^- A ^). 
k P 

C)=£ /C(S) ( A w--0 
k P 



(31) 



x a = 22 l T k P 
k P 



-^"ka-^pd 

x [sin A6 aa cos A(9 p k + cos A9 aa sin A# p k] 



(35) 



with A9 aa = Q\ a — 9^ a , and A kctQ is the surviving order pa- 
rameter. As expected, the spin-current changes sign depend- 
ing on whether it is the A k fT or A k xx order parameter that is 
present. 

For collinear magnetization = 0), an ordinary Josephson 
effect occurs with the superconducting phase difference as the 
driving force. Interestingly, one is able to tune both the spin- 
and charge-current to zero in the A 1 -phase when rriL || — mp> 
(0 = 7r). It follows from Eq. d34l ) that the spin- and charge- 
current only differ by a constant pre-factor 



tp/ -'tp,. 



-cue, 



a 



±1. 



(36) 
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It is then reasonable to draw the conclusion that we are dealing 
with a completely spin-polarized current such that both 
and J t p z must vanish simultaneously at i3 — n. 

Another result that can be extracted from Eqs. (1321 and (1341 
is a persistent non-zero DC spin-Josephson current even if the 
magnetizations on each side of the junction are of equal mag- 
nitude and collinear ($ = 0). This is quite different from 
the spin-Josephson effect recently considered in ferromag- 
netic metal junctions 24 . In that case, a twist in the magnetiza- 
tion across the junction is required to drive the spin-Josephson 
effect. 

Note that in the common approximation T kp = T, i.e. the 
tunneling probability is independent of the electron magnitude 
and direction of electron momentum, the two-particle current 
predicted above is identically equal to zero. Of course, such a 
crude approximation does not correspond to the correct phys- 
ical picture (see e.g. Ref. 51), and in general one cannot ne- 
glect the directional dependence of the tunneling matrix el- 
ement. This demonstrates that we are dealing with a more 
subtle effect than what could be unveiled when applying the 
approximation of a constant tunneling matrix element. 

An interesting situation arises in the case of zero externally 
applied voltage and identical superconductors on each side of 
the junction with SC phase differences A9 aa = 0. In this 
case, we find that Ij = while 

II = -2^|T kp | 2 sm 2 (tf/2)|A kn A pi jF kpU 
k P 

x sin(^ - 6f)/(E k1 E pl ). (37) 

when eV = 0, A8 acr = 0. Thus, we have found a dissi- 
pationless spin-current in the two-particle channel without an 
externally applied voltage and without a SC phase difference. 
This effect is present as long as is not or ir, correspond- 
ing to parallel or anti-parallel magnetization on each side of 
the junction. It is seen from Eq. (137) 1 that the spin-current is 
driven by an interband phase difference on each side of the 
junction. A necessary condition for this effect to occur is that 
no inter-band Josephson coupling is present, i.e. electrons in 
the two energy-bands and _E k j do not communicate with 
each other. To understand why a Josephson coupling would 



destroy the above effect, consider the free energy density for 
a p-wave FMSC first proposed in Ref. given by 

^JT'-Ajcos^ -0 U ) (38) 

in the presence of a Josephson coupling. In Eq. d38l . Aj deter- 
mines the strength of the interaction while T' contains the SC 
and FM contribution to the free energy density in addition to 
the coupling terms between the SC and FM order parameters. 
Consequently, the phase difference 6^ — On is locked to or 
7r in order to minimize T, depending on sgn(Aj). Considering 
Eq. d37l . we see that if z = in this case, since the argu- 
ment of the last sine is zero. Mechanisms that would induce 
a Josephson coupling include magnetic impurities causing in- 
elastic spin-flip scattering between the energy-bands and spin- 
orbit coupling. Recently, the authors of Ref. l50l proposed that 
p-wave SC arising out of a FM metal state could be explained 
by the Berry curvature field that is present in ferromagnets 
with spin-orbit coupling. It is clear that in the case where spin- 
orbit coupling is included in the problem, spin-flip scattering 
processes occur between the energy bands such that the f and 
I spins can not be considered as two independent species any 
more. The SC phases will then be locked to each other with 
a relative phase of or w. However, note that in the general 
case, Eq. d32l produces a non-zero charge- and spin-current 
even if the spin-up and spin-down phases are locked to each 
other. 



E. Single-particle currents 

In the single-particle channel, we find that the charge and 
spin-currents read 

4,« = E a <^«"W>^ (39) 

a 

as seen from Eq. ( f29b . From Eq. ( lAlb . we then extract the 
proper expectation value, which is found to be 



{N aa (t)), P = 4tt^[1 + aacos(tf)]|T kp 

a | ^ka"^j 

kpa" 



2 

pa 



[/(£ka) - f{E pa )] (6{-eV + E ka - E v „) 



ikaa^pcrcj | 



(fka + E ka ) 2 (£ pa +Epa 



; S(-eV - E pa + E ka ) 



[1 - /(^ka) - /(£, 



pa) 



Ak '" ! '~ ;5(-eV - E ka - E pa ) - - } Ap °f 5(-eV + E ka + E pa ) 



(Cka +E ka ) 2 



(Cpo- + Epg) 



(40) 



The currents in Eq. d39b are thus seen to require an applied this is because the Cooper pairs need to be split up in order 
voltage in order to flow in the tunneling junction. Clearly, 
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for a single-particle current to exist, such that both spin- and 
charge-currents vanish at eV = 0. 

In Ref. [iU, the presence of a persistent spin-current in the 
single-particle channel for FM/FM junctions with a twist in 
magnetization across the junction was predicted. For consis- 
tency, our results must confirm this prediction for the single- 
particle current in the limit where SC is lost, i.e. Ak CTCT — > 0. 
Note that the z-direction in Ref. [24] corresponds to a vector 
in our local :rj/-plane since the present quantization axis lies 
parallel with the magnetization direction. Upon calculating 
the x- and y-components of the single-particle spin-current 
for our system in the limit where SC is lost, i.e. A^ aa — > 0, 
a persistent spin Josephson-like current proportional to sin($) 
is identified. More precisely, 

kp after 

><M<W^(- e ^)} (41) 

when A q(T(J = (see Appendix for details). In agreement 
with Ref. [2J, the component of the spin-current parallel to 
niL x rriR is seen to vanish for •& — {0, 7r} at eV = 0. 



III. FERROMAGNETS WITH SPIN-ORBIT COUPLING 

A. Coexistence of ferromagnetism and spin-orbit coupling 

In a system where time-reversal and spatial inversion sym- 
metry are simultaneously broken, it is clear that spins are 
heavily affected by these properties. There is currently 
much focus on ferromagnetic semiconductors where spin- 
orbit coupling plays a crucial role with regard to transport 
properties 21 ' 22 . In fact, there has in recent years been much 
progress in the semiconductor research community where the 
spin-Hall effect in particular has received much attention 23 . 
With the discovery 5 - 2 - of hole-mediated ferromagnetic order 
in (In,Mn)As, extensive research on III-V host materials 
was triggered. Moreover, it is clear that properties such 
as ferromagnetic transition temperatures in excess of 100 
KM and long spin-coherence times 54 in GaAs have strongly 
contributed to opening up a vista plethora for information 
processing and storage technologies in these new magnetic 
mediums 55 . 

Generally, spin-orbit coupling (SOC) can be roughly di- 
vided into two categories - intrinsic and extrinsic. Intrinsic 
SOC is found in materials with a non-centrosymmetric crystal 
symmetry, i.e. where inversion symmetry is broken, whereas 
extrinsic SOC is due to asymmetries caused by impurities, lo- 
cal confinements of electrons or externally applied electrical 
fields. 

In the present paper, we investigate the tunneling current of 
spin between two ferromagnetic metals with spin-orbit cou- 
pling induced by an external electric field. This way, we will 
have two externally controllable parameters; the magnetiza- 
tion m and the electrical field E. The case of tunneling be- 
tween two noncentrosymmetric superconductors with signif- 



icant spin-orbit coupling, but no ferromagnetism, has previ- 
ously been considered in Ref. H^. 

B. The Hamiltonian 

Our system consists of two Heisenberg ferromagnets with 
substantial spin-orbit coupling, separated by a thin insulating 
barrier which is assumed to be spin-inactive. This is shown in 
Fig. |4] We now operate with only one quantization axis, such 
that a proper tunneling Hamiltonian for this purpose is 

Hj = ^(T kp cLd pCT + h.c), (42) 

kpo" 

where {c ko .,Ck<r} and {c2 kcr , d^} are creation and annihila- 
tion operators for an electron with momentum k and spin a 
on the right and left side of the junction, respectively, while 
Tk p is the spin-independent tunneling matrix element. In k- 



m L Thin, insulating barrier 
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.-' FM with S-0 coupling 
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FIG. 4: Our model consisting of two ferromagnetic metals with spin- 
orbit coupling separated by a thin insulating barrier. The magnetiza- 
tion m and electrical field E are allowed to point in any direction so 
that our results are generally valid, while special cases such as planar 
magnetization etc. are easily obtained by applying the proper limits 
to the general expressions. 

space, the Hamiltonian describing the ferromagnetism reads 
#fm = J2 £ kcL c k. - »7(k)S k • S_ k (43) 

kcr k 

in which £k is the kinetic energy of the electrons, J is the 
ferromagnetic coupling constant, N is the number of parti- 
cles in the system, while S k = (1/2) ^ Q(3 ^La^ap^ is the 
spin operator. As we later adopt the mean-field approxima- 
tion, m = (m x , m y , m z ) will denote the magnetization of the 
system. 

The spin-orbit interactions are accounted for by a Rashba 
Hamiltonian 

ffs-o = - x k) • <% k , (44) 

k 

where <^k = [ckf, Cki] T , E = — VU is the electrical field felt 
by the electrons and <x = (di, 62, 03) in which <jj are Pauli 
matrices, while the parameter £ is material-dependent. From 
now on, the notation £(E x k) = Bk = {Bk >x , B^ y, i?k,z) 
will be used. In general, the electromagnetic potential V con- 
sists of two parts V ml and V ext (see e.g. Ref. |23| for a detailed 
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discussion of the spin-orbit Hamiltonian). The crystal poten- 
tial of the material is represented by V[ nt , and only gives rise 
to a spin-orbit coupling if inversion symmetry is broken in the 
crystal structure. Asymmetries such as impurities and local 
confinements of electrons are included in V ext , as well as any 
external electrical field. Note that any lack of crystal inver- 
sion symmetry results in a so-called Dresselhaus term in the 
Hamiltonian, which is present in the absence of any impurities 
and confinement potentials. In the following, we focus on the 
spin-orbit coupling resulting from V ext , thus considering any 
symmetry-breaking electrical field that arises from charged 
impurities or which is applied externally. In the case where 
the crystal structure does not respect inversion symmetry, a 
Dresselhaus term 57 can be easily included in the Hamiltonian 
by performing the substitution 



(E x k) • & -> [(E x k) + D(k)] • &, 



(45) 



where V(k) = -D(-k). 

We now proceed to calculate the spin-current that is gener- 
ated across the junction as a result of tunneling. Note that in 
our model, the magnetization vector and electrical field are al- 
lowed to point in arbitrary directions. In this way, the obtained 
result for the spin-current will be generally valid and special 
cases, e.g. thin films, are easily obtained by taking the appro- 
priate limits in the final result. It should be mentioned that the 
effective magnetic field from the spin-orbit interactions might 
influence the direction of the magnetization in the ferromag- 
net. This is, however, not the main focus of our work, and 
we leave this question open for study. Our emphasis in the 
present paper concerns the derivation of general results onto 
which specific restrictions may be applied as they seem ap- 
propriate. 

In the mean-field approximation, the Hamiltonian for the 
right side of the junction can be written as H = i/pM + H$_o, 
which in a compact form yields 



H - 



k 



4 + B k H 



-Cr + B K - 

£ ki 



<Pk, (46) 



where e k(T = £k — f (Cz,R — B k . z ) and Ho is an irrelevant 
constant. The FM order parameters are Cr = 2 Ji](0)(itir^ x — 
im RiV ) and (z.r = 2J7?(0)m R , z and B k ,± = B k>x ± i-Bk, y - 
For convenience, we from now on write C = |C|e'^ and 
Bk,± = \Bk t ±\e^ 1Xk . The Hamiltonian for the left side of 
the junction is obtained from Eq. d46b simply by the doing the 
replacements k — > p and R — » L. 



C. Tunneling formalism 

In order to obtain the expressions for the spin- and charge- 
tunneling currents, it is necessary to calculate the Green func- 
tions. These are given by the matrix 

= (-iu> n l + A k )-\ (47) 

where _4 k is the matrix in Eq. (l46l i. Explicitly, we have that 



(48) 



Above, u> n = 2(n + l)ir/ (3,n = 0, 1, 2 ... is the fermionic 
Matsubara frequency and (3 denotes inverse temperature. In- 
troducing 

X k (iu) n ) = (e kT - iw„)(ekj. -iu n ) ~ |Cr - -Bk,-| 2 , (49) 

the normal and anomalous Green functions are 

G k ff (iw„) = (e k _ CT - iw„)/X k (icj n ), 

F k iT (i^„) = F^'\ju n ) = (Cr - S k) _)/X k (iw„). (50) 

The expression for I s (t) is established by first consider- 
ing the generalized number operator N a p = ^ k c tQ c k ^. 
This operator changes with time due to tunneling according to 
N a p = i[-Ht, Nap], which in the interaction picture represen- 
tation becomes N afj (t) = -iI] k p(7 kp c ka dp/3e lteV ' -h.c). 
The voltage drop across the junction is given by the difference 
in chemical potential on each side, i.e. eV = /jr — /iL- In the 
linear response regime, the spin-current across the junction is 



i s (f) 



(51) 



a/3 



where the expectation value of the time derivative of the trans- 
port operator is calculated by means of the Kubo formula Eq. 
(l30l >. Details will be given in Sec. IA2I 



D. Single-particle currents 

At eV — 0, it is readily seen from the discussion in 
Sec. IA2I that the charge-current vanishes. Consider now 
the z-component of the spin-current in particular, which can 
be written as if = 3m{$(- eV)}. The Matsubara func- 
tion <£>(— eV) is found by performing analytical continuation 
— » — eV + i0 + on $(£;„), where 

$(£>„) = - I>(Gr(^m)G7(i^ m -i£,) 

+F- a >°(uj m )F°>-°(iu m - £;„)) . 

(52) 

Here, uj v — 2vn/ f3, v = 0, 1, 2 ... is the bosonic Matsubara 
frequency. Inserting the Green functions from Eq. (l50l into 
Eq. (l52l . one finds that a persistent spin-current is established 
across the tunneling junction. For zero applied voltage, we 
obtain 



E 



^kpl^kp 



kp 



27k7 P 



CrCl| sinA</> + \B k -B p -\ sin Axk P 
S k _Cl| sin(x k - L ) - l-Bp _Cr| sm{c/) R - x P ) , 
n(e k + a7 k ) - n(e p + /?7 P ) 



a = ± 

|3=± 



(e k + «7 k ) - (e p + /?7 P ) 



(53a) 
(53b) 
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In Eqs. ||53), A^kp = Xk - X P , A<p = L , while 

7k = (C,,R-S k , 2 ) 2 + |C R -B k ,-| 2 (54) 

and n(e) denotes the Fermi distribution. In the above expres- 
sions, we have implicitly associated the right side R with the 
momentum label k and L with p for more concise notation, 
such that e.g. B^ z = -B k Defining Q — 2Jrj(0)m,i, we see 
that Eq. d54"l i can be written as 

7k = |C R -B k |. (55) 

The spin-current described in Eq. d53l i can be controlled by 
adjusting the relative orientation of the magnetization vectors 
on each side of the junction, i.e. A</>, and also responds to 
a change in direction of the applied electric fields. The pres- 
ence of an external magnetic field H; would control the ori- 
entation of the internal magnetization irij . Alternatively, one 
may also use exchange biasing to an anti-ferromagnet in order 
to lock the magnetization direction. Consequently, the spin- 
current can be manipulated by the external control parameters 
{Hi, Ei} in a well-defined manner. This observation is highly 
suggestive in terms of novel nanotechnological devices. 

We stress that Eq. ([53} is non-zero in the general case, since 
7k 7^ — 7-k and x_k = Xk + tt- Moreover, Eq. d53l is 
valid for any orientation of both m and E on each side of the 
junction, and a number of interesting special cases can now 
easily be considered simply by applying the appropriate limits 
to this general expression. 



such that the summation over a and (3 yields zero. Thus, no 
spin-current will exist at eV — over a tunneling barrier sep- 
arating two systems with spin-orbit coupling alone. In the 
general case where both ferromagnetism and spin-orbit cou- 
pling are present, the density of states at, say, Fermi level are 
different, leading to a persistent spin-current across the junc- 
tion due to the difference between (E) and (E). 



W W M 




o o o 



FIG. 5: Schematic illustration of the energy-bands for (a) a system 
with spin-orbit coupling, (b) a system with ferromagnetic ordering, 
and (c) a system exhibiting both of the aforementioned properties. 
The dotted line corresponds to quasi-particles with a =\ r , while 
the full drawn line designates a =f. Since the density of states 
N" (Eka-) is proportional to (dEka/dk)" 1 , we see that a difference 
between (Ek&) and (Ek&) is zero at Fermi level in (a), while 
the density of states differ for the j- and .[.-populations in (b) and (c). 
Thus, a persistent spin-current will only occur for tunneling between 
systems corresponding to (b) and (c). 



E. Special limits 

Consider first the limit where ferromagnetism is absent, 
such that the tunneling occurs between two bulk materials 
with spin-orbit coupling. Applying m — > to Eq. (|53l , it 
is readily seen that the spin-current vanishes for any orienta- 
tion of the electrical fields. Intuitively, one can understand 
this by considering the band structure of the quasi-particles 
with energy E^ a = £ k + C7k and the corresponding density 
of states N(Ekcr) when only spin-orbit coupling is present, as 
shown in Fig. [5] Since the density of states is equal for f and 
| spins 69 , one type of spin is not preferred compared to the 
other with regard to tunneling, resulting in a net spin-current 
of zero. Formally, the vanishing of the spin-current can be un- 
derstood by replacing the momentum summation with integra- 
tion over energy, i.e. £ kp -> J J dE R dE L N R {E R )N L {E L ). 
When m — > 0, Eq. ([53j dictates that 



r 



E 



a(3 



dE K!a dE h>0 Ng{E K , a )Nl(E ht0 ) 



(56) 



Since the density of states for the f- and | -populations are 
equal in the individual subsystems, i.e. (E) = N^-(E) = 
N(E), the integrand of Eq. d56l l becomes spin-independent 



We now consider a special case where the bulk structures 
indicated in Fig. |4] are reduced to two thin-film ferromag- 
nets in the presence of electrical fields that are perpendicular 
to each other, say E L = (E h ,0, 0) and E R = (0, Er, 0), as 
shown in Fig. |6ja) and (b). In this case, we have chosen 
an in-plane magnetization for each of the thin-films. Solving 
specifically for Fig. |6fa), it is seen that rxiL = (0, m^j,, tol,z) 
and rriR = (mR )X , 0, tor iZ ). Furthermore, assume that the 
electrons are restricted from moving in the "thin" dimension, 
i.e. p = (Q,p v ,Pz) and k = (fc x , 0, k z ). In this case, Eq. (l53l l 
reduces to the form 

J z = ^osgn(m L ,j/) + 7 i,k P sgn(p z ), (57) 
k P 

where the constants above are 

|T k p| 2 J kp (|CRCL|-£R|fc z CL|) 



'o = £ 



k P 



2|CR + B k ||C L + Bp| 



A.kp — 



|Tkp| 2 Jk P £ L (£R|fc z p z | - | Pz Cr|) 
2|C R + Bk||C L + Bp| 



(58) 



with 



Cl = 2j7y(0)(0,m L ,y,m L , 2 ) 
C R = 2Jj7(0)(mR, x ,0,mR, z ) 



B p = £ h E h (0,-p z ,p y ) 
B* = Z R E K (k v ,a,-k x ) 



(59) 
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such that ii.kp ^ -Z^-k.-p- Likewise, for the setup sketched 
in Fig. |6|b), one obtains 



|7k p | 2 Jkp 



kp 



2|Cr + E R (k y - ifc K )| 2 |CL + E h ( Py - i Px )\ 



x [KrClI sin + E R E h (k 2 y + k 2 x )(p 
-E R \C L \(kl + kl)sm( Xk -M 
~ El|CrI(p£ + ^)sin(0 R - x P ) 

where x q obeys 



2 ' Px)|sinAxk P 



tanxq 



<7.r 



, q = k,p. 



(60) 



(61) 



From these observations, we can draw the following con- 
clusions: whereas the spin-current is zero for the system in 
Fig. |6|a) and (b) if only spin-orbit coupling is considered, it 
is non-zero when only ferromagnetism is taken into account. 
However, in the general case where both ferromagnetism and 
spin-orbit coupling are included, an additional term in the 
spin-current is induced compared to the pure ferromagnetic 
case. Accordingly, there is an interplay between the magnetic 
order and the Rashba-interaction that produces a spin-current 
which is more than just the sum of the individual contribu- 
tions. 
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FIG. 6: Tunneling between planar ferromagnets in the presence of 
externally applied electrical fields El and Er that destroy inversion 
symmetry and induce a spin-orbit coupling. 



IV. DISCUSSION OF RESULTS 

Having presented the general results for tunneling currents 
between systems with multiple broken symmetries in the pre- 
ceeding sections, we now focus on detection and experimental 
issues concerning verification of our predictions. 

Consider first the system consisting of two ferromagnetic 
superconductors separated by a thin, insulating barrier. It is 
well-known that for tunneling currents flowing between two 
s-wave SC in the presence of a magnetic field that is perpen- 
dicular to the tunneling direction, the resulting flux thread- 
ing the junction leads to a Fraunhofer-like variation in the DC 
Josephson effect, given by a multiplicative factor 



sin(7T$/$ ) 
(tt$/$ ) 



(62) 



in the critical current. Here, $o = nti/e is the elementary flux 
quantum, and $ is the total flux threading the junction due 
to a magnetic field. Consequently, the presence of magnetic 
flux in the tunneling junction of two s-wave SC threatens to 
nullify the total Josephson current. In the present case of two 
p-wave FMSC, this is not an issue since we have assumed uni- 
form coexistence of the SC and FM order parameters which is 
plausible for a weak intrinsic magnetization. The effect of an 
external magnetic field H would then simply be to rotate the 
internal magnetization as dictated by the term H • m in the 
free energy T (see e.g. Ref.lla). Thus, there is no diffraction 
pattern present for the tunneling-currents between two non- 
unitary ESP FMSC, regardless of how the internal magnetiza- 
tion is oriented. Since the motion of the Cooper-pairs is also 
restricted by the thin-film structure, there is no orbital effect 
from such a magnetization. 

Note that the interplay between ferromagnetism and su- 
perconductivity is manifest in the charge- as well as spin- 
currents, the former being readily measurable. Detection 
of the induced spin-currents would be challenging, although 
recent studies suggest feasible methods of measuring such 
quantities 58 . We comment more on this later in this sec- 
tion. First, we adress the issue of how boundary effects af- 
fect the order parameters. Studies 40 ' 41,42 have shown that in- 
terfaces/surfaces may have a pair-breaking effect on uncon- 
ventional SC order parameters. This is highly relevant in 
tunneling junction experiments as in the present case. The 
suppression of the order parameter is caused by a forma- 
tion of so-called midgap surface states (also known as zero- 
energy states) 39 which occurs for certain orientations of the 
k-dependent SC gaps that satisfy a resonance condition. Note 
that this is not the case for conventional s-wave superconduc- 
tors since the gap is isotropic in that case. This pair-breaking 
surface effect was studied specifically for p-wave order pa- 
rameters in Refs. and it was found that the component 
of the order parameter that experiences a sign change under 
the transformation k± — > — k±, where k± is the component of 
momentum perpendicular to the tunneling junction, was sup- 
pressed in the vicinity of the junction. By vicinity of the junc- 
tion, we here mean a distance comparable to the coherence 
length, typically of order 1-10 nm. Thus, depending on the 
explicit form of the superconducting gaps in the FMSC, these 
could be subject to a reduction close to the junction, which in 
turn would reduce the magnitude of the Josephson effect we 
predict. Nevertheless, the latter is nonvanishing in the general 
case. 

Since the critical Josephson currents depend on the rela- 
tive magnetization orientation, one is able to tune these cur- 
rents in a well-defined manner by varying This can be 
done by applying an external magnetic field in the plane of 
the FMSC. In the presence of a rotating magnetic moment 
on either side of the junction, the Josephson currents will 
thus vary according to Eq. d32l . which may be cast into 
the form if = Iq + / m cos(#). Depending on the relative 
magnitudes of Iq and I m , the sign of the critical current may 
change. Note that such a variation of the magnetization vec- 
tors must take place in an adiabatic manner so that the sys- 
tems can be considered to be in, or near, equilibrium at all 
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times. Our predictions can thus be verified by measuring the 
critical current at eV = for different angles d and com- 
pare the results with our theory. Recently, it has been re- 
ported that a spin-triplet supercurrent, induced by Josephson 
tunneling between two s-wave superconductors across a fer- 
romagnetic metallic contact, can be controlled by varying the 
magnetization of the ferromagnetic contact 59 . Moreover, con- 
cerning the spin- Josephson current we propose, detection of 
induced spin-currents are challenging, although recent stud- 
ies suggest feasible methods of measuring such quantities 58 . 
Observation of macroscopic spin-currents in superconductors 
may also be possible via angle resolved photo-emission exper- 
iments with circularly polarized photons 60 , or in spin-resolved 
neutron scattering experiments 32 . 

We reemphasize that the above ideas should be experimen- 
tally realizable by e.g. utilizing various geometries in order to 
vary the demagnetization fields. Alternatively, one may use 
exchange biasing to an anti-ferromagnet. Such techniques of 
achieving non-collinearity are routinely used in ferromagnet- 
normal metal structures—. 

With regard to the predicted DC spin-current in for a system 
consisting of two ferromagnetic metals with spin-orbit cou- 
pling, we here suggest how this effect could be probed for in 
an experimental setup. For instance, the authors of Ref. |62j 
propose a spin-mechanical device which exploits nanome- 
chanical torque for detection and control of a spin-current. 
Similarly, a setup coupling the electron spin to the mechani- 
cal motion of a nanomechanical system is proposed in Ref.HH 
The latter method employs the strain-induced spin-orbit inter- 
action of electrons in a narrow gap semiconductor. In Ref.l63l 
it was demonstrated that a steady-state magnetic-moment cur- 
rent, i.e. spin-current, will induce a static electric field. This 
fact may be suggestive in terms of detection 64 ' 65 , and could be 
useful to observe the novel effects predicted in this paper. 



V. SUMMARY 

In summary, we have considered supercurrents of spin and 
charge that exist in FMSC/FMSC and FMSO/FMSO tunnel- 
ing junctions. In the former case, we have found an interplay 
between the relative magnetization orientation on each side 
of the junction and the SC phase difference when considering 
tunneling between two non-unitary ESP FMSC with coexist- 
ing and uniform FM and SC order. This interplay is present in 
the Josephson channel, offering the opportunity to tune dissi- 
pationless currents of spin and charge in a well-defined man- 
ner by adjusting the relative magnetization orientation on each 
side of the junction. As a special case, we considered the 
case where the SC phase difference is zero, and found that a 
dissipationless spin-current without charge-current would be 
established across the junction. Suggestions concerning the 
detection of the effects we predict have been made. 

Moreover, we have derived an expression for a dissipa- 
tionless spin-current that arises in the junction between two 
Heisenberg ferromagnets with spin-orbit coupling. We have 
shown that the spin-current is driven by terms originating from 
both the ferromagnetic phase difference, in agreement with 



the result of Ref. 24, and the presence of spin-orbit coupling 
itself. In addition, it was found that the simultaneous breaking 
of time-reversal and inversion symmetry fosters an interplay 
between ferromagnetism and spin-orbit coupling in the spin- 
current. Availing oneself of external magnetic and electric 
fields, our expressions show that the spin-current can be tuned 
in a well-defined manner. These results are of significance 
in the field of spintronics in terms of quantum transport, and 
offer insight into how the spin-current behaves for nanostruc- 
tures exhibiting both ferromagnetism and spin-orbit coupling. 
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APPENDIX A: DETAILS OF MATSUBARA FORMALISM 
1. Ferromagnetic superconductors 

Inserting Eq. d27l ) into Eq. d30l ), one finds that 

(N a p(t)) = (iV a/ j(t)) sp + (N a p(t)) lp 



dt 



([M a p(t),Mi(t')])i 



-ieV(t-t') 



{[Ml a (t),M{t')])^ v ^ 

([M af3 (t),M(t')])e- leV ^ 

([MUt),M^t')]}e ie nt+t') 



(Al) 



where the two first terms in Eq. ( IAU contribute to the single- 
particle current while the two last terms constitute the Joseph- 
son current. Above, we defined 

M a(j {t) = ^^WTkpcLWVW 

kptr 

M(t)= ^iv 2 V) T k P cLWV'W- (A2) 

By observing that <7 Q/3 = (cr^)*, we can combine Eqs. (130b - 
(TA2l to yield 

& af} (Naf3{t)U = 23m{er Q/3 $ Q/3 , sp (- e y)} 
*a/s<^»(*))tP = ^m{& a ^ aP j(eV)e- 2letv } (A3) 

where the Matsubara functions are obtained by performing an- 
alytical continuation according to 

$a/3,sp(-eV) = _ lim $ o /3,sp(l£5i0 
uo w — > — eV+iO+ 

®a/3,](eV) = _ lim $ a/3 ,u,(iav), (A4) 
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In Eq. ( IA4I ). ui„ — 2wu/j3, v — 1, 2, 3 . . . is the bosonic 
Matsubara frequency and 



Jo 



We then Fourier-transform Eq. (|A9t into 



G (k,r) = l^c- lu mT G a (p, ] 



(A10) 



kp 



Tk P T k % p , (T{4 a (r) d pCT (t)^ , CTi (0)c k , CT2 (0) }} , 
^(^) = - / dre 1 ^ J] ^ /2 V)^Wx 

•'O 

kper 
k'p'o-i ct 2 

T kp T kV (f{4 Q (r)V«4v 2 (0)V.i(0)})- 

(A5) 

Here, T denotes the time-ordering operator, and (3 = 
l/fc B T is the inverse temperature. Only k' = (— )k, p' = 
(— )p contributes in the single-particle (two-particle) chan- 
nel, while the diagonalized basis <j5 k(T dictates that only ctq, = 
a, u\ = a contributes in the spin summation. Making use of 
the relation k(T = (p^Uua, Eq. (1A51 > becomes 



where uj m = (2m + l)n//3, m — 1, 2, 3 ... is a fermionic 
Matsubara frequency. The frequency summation over m is 
evaluated by contour integration as in e.g. Ref. 66 to yield the 
result 



1 V^a, /, n,-,/ ~ , s f(Eka) — f{Epcr) 

- > G a (k.,w m )G a (p,vjj u +w m ) = — — — 

-^Ga^LO^Ga^^Uy-lUJrn) = ~ _ ^ — p 







1U) U — Efra — Epv 

(All) 



k P <T 

x (n[Un^L(T) + ^2ka7-ka(T)" 
*7llka7ka(0) + J7i 2kQ 7i kQ (0)] }} 

Unp*lU°) + tfiWr-k^O)! }) (A6) 



where f(E) = l- f(-E) = 1/(1 + c^) is the Fermi distri- 
bution. It is then a matter of straight-forward calculations to 
obtain the result 



kpo" 



2 

pa 



(£ka + -Eka)(£p<T + E pa ) kp<T 

, |Ap CTff | 2 A kp i(-eF) 



Jo . 



kpcra V 



lAkaapAiL^-ey) 



£ka + Ska 



(A12) 



kpcr 

„t 



X t/21ka7k CT (0) + L>22ka7i kc7 (0) }) 

x (t[?/ 2 * 1pct7pct (o) + U* 22pa i-pM 

X [f/llpa7p<r(T) + C>12p CT 7ip CT (T)]}> (A7) 

Since our diagonalized Hamiltonian has the form of a free- 
electron gas, i.e. 



1W(eF) = - V% 2) {d)V^ (tf)T kp T_k,- 



kpo" 



f-^ E A k PCTtt (^), (Ai3) 



4i? kQ i? p(T 



A=±l 
p=±l 



where A kpcrQ (eV r ) is obtained by performing analytical con- 



H F MSC = H + Y £; k<T 7 k CT 7k 



kp< 

tinuation \uj v — > + i0 + 



on 



(A8) 



ktr 



with _ffo = -Ho — (^kt + ^kj.)> the product of the new fermion 
operators <^ k(T in Eq. (IA6t yield unperturbed Green's func- 
tions according to 

G a (k, t-t') = (T{ c t (r')c ka (r)}) (A9) 



A[/(Ska)-/(A^ pCT )], 



VjJ v + pEk a — \E pa 



X,p = ±l. 



(A14) 

We also provide the details of the persistent spin- 
supercurrent for A aa = 0. Writing the Josephson current 
Eq. d32l out explicitly, one has that if = el + and if = — I~ 
where 
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k P 

+ sin 2 (i9/2) 



C0S 2 W 2)l^n|pH sinA%Fkpn 



-E-kT-E-pT 



l A kTT A plll 



7 U 



± sin 2 W2) |A ^^ pTTl sin(^ T - 0* )F kpiT 
i'ki-fe'pT 



±cos 2 0V2) 



l A k|jA pU | 



sinA0 u F kpU 



(A15) 



The first and fourth term above vanish when A8 aa = 0. By 
observing that -Fkptl = -Ppkj.T> we 816 tnen able to re-write 
Eq. dA15t as 

^=£|Tk P | 2 sin 2 W2)^l^^ pn 



x 



[sm(0j;-0f)±sin(0^-0^ 



= e £|T kp | 2 sin 2 (tf/2) 
k P 



|A WT A 



kT^pl 



^kf-Epi 

[sin^^T^i -Bf ±0\)/2 
cos((^±6'f-^T6 , T)/2 



F) 



kpTJ 



(A16) 



It is clear that the argument of the sine gives for the upper 



J 



sign, such that if = 0. But for the lower sign, the argument 
of the cosine is equal to 0, such that Eq. ( |37| i is obtained. 



2. Ferromagnets with spin-orbit coupling 

The spin-current across the junction can be written as 

I s = 3fm{*(-eV)}, 



<&(— eV) = lim ^(iuv), 



where we have defined the Matsubara function 



(A17) 



*(iS5„) = £ |T kp | 2 <x Q/3 / dre^ 

kpa/3<r ^° 

x (T{c kCT (0)cL (r) }) (T{d p/3 (r)d pff (0) }) . 

(A18) 



In Eq. JA18t . we defined the time-ordering operator T while 
/3 in the upper integration limit is inverse temperature and 
ui v = 2nir/j3, n = 0, 1, 2, ... is a bosonic Matsubara fre- 
quency. From the definition of the spin-generalized Green's 
function 

Gf (r t>) = -(T{cUt)cI p (t>)}), (A19) 
Eq. JA1 8I > can be written out explicitely to yield 



*(i<3„) = i £ |T Up | 



kp,m 



■ n yG^ (icj m )G P T (jw m - + G k T (iuj m )G P l (iuj m -H2 v ) y j+& n (G^ (iw m )G p T (l 
+ G^ (iuj m )G p l (iuj m - iu^)) + o- iT ^Gj c i (iw m )G P T (iw m - iC) + Gl L {iuj m )G p L {iLu m - i2„)J 
+ <r u (G k x (iw m )G p T (iw m - iw„) + G k x (jw m )G^ (jw m - £>„)) . (A20) 



We made use of the Fourier-transformations 

Gf(icj m )= [ dre^Gf (r), 
Jo 



(A21) 



in writing down Eq. (1A201 >. where oj m — 2(m + l)ir/(3, m — 
0, 1, 2, ... is a fermionic Matsubara frequency. Having written 
down the full expression for the Matsubara function in Eq. 
< IA20b . one can now easily distinguish between components of 
the spin-current. For instance, only a aa will contribute to the 
z-component of I s , and the corresponding terms can be read 
out from Eq. ( 1A20I ). From the present Green functions in Eq. 



r 



( T50l ), it is obvious that three types of frequency summations 
must be performed, namely 



Jkp,r — g 



[(e kT - iw m )(e k j. - iw m ) - vl] 
1 

[(e pT - tw m + iS i/ )(e P 4 - tw m + v2 v ) - y p ] 

(A22) 
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with r is an integer. Performing the summation over m using with the definition ^ ka = £k + o.y\s.- Separating the gen- 
residue calculus, one finds that eral expression Eq. dA2()j l into its spatial components * = 

(Q x , the components of the spin-current can be ex- 

tracted according to If = Sm{$i(-ey)}, i — x,y,z. 
Note that the charge-current in this model, which vanishes 
for eV = 0, is obtained by the performing the replacement 
(A23) a a j3 — > l a p, where 1 is the 2x2 unit matrix. We find that 



"kp,r 



a(3 



y 

/3=± 



k P 



i^kpi 2 

47k7 P 



Jk P ,o(ek|(CL - Bp,-) + (e pT + - B k ,_) + (e P l + i^XCr - B k ,+) + e kT (4 - Bp,-) 

- Jk P ,i((C L - B P ,-) + (Cr - B k ,_) + (4 - B k ,+) + - B Pl+ ) 



IT, 



kp| 



kp 



47k7 P 



IT, 



kp| 



kp 



47k7 P 



Jk P ,of - £ki(CL - Bp,-) - (e P t +i^)(Cr - B k) _) + (e pi + £i v )((l - B K+ )+e^{(l - B p ,_) 
- Jk P ,i ( - (Cl - B p ,_) - (Cr - B k ,_) + (4 - B k .+) + (4 - B Pr 

</k P ,o(£ki(£pi + - £kT(e P T + iWf) + (Cr - B k ,_)(4 - B p ,+) - (C R - B kj+ )(C L - B p 



+ '/ kp .i (tTcT - £ki + £ P T - £p;J 



(A24) 
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